Squeezing mode in nuclear collisions 
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The time dependent Schroedinger equation is solved analytically for a simplified model of moving 
infinite walls. A new knock-out mode is described which might occur during heavy ion collisions. The 
outer shell-nucleons are ionised due to the increase of level energy when two nuclei are approaching 
fast enough. This is analogous to the Mott effect but in contrast occurs only if the reaction time 
is short enough that no common ionisation threshold in the compound system is established. To 
demonstrate this pure nonequilibrium effect a simulation of realistic heavy ion collision by a nonlocal 
Boltzmann equation is performed. 
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When a system with bound states is exposed to a com- 
pression beyond certain values the bound states break off 
and decompose into their constituents. This pressure ion- 
isation known as the Mott transition is well established in 
different fields of physics Alternatively in a many 

body system the density could be increased. The theo- 
retical treatments agree in that the ionisation threshold 
is lowered faster than the binding energy with increas- 
ing density which leads to a cross-over and ionisation at 
the Mott density. These treatments rely on the fact that 
one has a certain degree of equilibration in the system, at 
least that the system has one unique ionisation threshold. 

This situation is however somewhat different when pro- 
cesses occur far from equilibrium. Then there might be 
not enough time to establish a common threshold in the 
system. In particular if two nuclei approach each other 
in a heavy ion collision it takes a certain time before a 
compound system is established or decomposition hap- 
pens at higher energies. One can easily imagine that 
there will be no common ionisation threshold for the 2 
nuclei at the early stage of reaction. Instead we will 
show that this leads to a new escaping mode by squeezing 
states which should lead to the nonequilibrium emission 
of particles. The principle phenomena has already been 
investigated in the past as light nonequilibrium particle 
emission Here we will show that with the help of an 
exactly solvable model of time dependent Schroedinger 
equation a new effect arise not described so far. The dif- 
ferent features are transversal angular distribution and 
a lower bound of projectile energy that this squeezing 
mode to happen. There arc experimental signals for dy- 
namical particle emission |||7|. In contrast to diabatic 
emission of particles Q limited to beam energies much 
below the Fermi energy we consider here the case of faster 
processes around the Fermi energy. After presenting the 
exactly solvable model we will confirm this mode by real- 
istic simulations which will show a transversal distribu- 
tion and low energy of emitted particles while diabatic 
emitted particles are longitudinal peaked. 

Another intuitive picture is the following. The Pauli 



principle will forbid over-occupation of states, which 
should result in a Fermi gap which is closed during dis- 
sipation. Before this quasi - equilibration happens one 
has essentially the situation that an outer nucleon feels 
a rapidly increasing force due to the other nuclei and 
the Pauli - forbidden areas in phase space. If the speed 
of nuclei is high there is no many - body equilibration 
or dissipation but rather a shrinking of phase space for 
the outer nucleons during the time after first touching 
of the two nuclei. Therefore it is reasonable to assume 
a picture where the boundary of the outer nucleons is 
shrinking with time. We will solve such a time depen- 
dent Schroedinger equation to show that indeed the level 
energy of the outer nucleons goes up and eventually leads 
to ionisation. Since this is opposite to the Mott transi- 
tion described above where the threshold decreases we 
call this mode squeezing mode in the following. There 
are similarities to the elevator resonant activation mode 
where a time dependent potential inside the wall cre- 
ates resonant levels which trap the particles and lift them 
above the barrier. Here we merely squeeze the wall. 

Let us assume that the outer shell nucleons are bound 
states which can be parameterised by a simple one- 
dimensional infinite wall model, i.e. a free particle inside 
an infinite wall of distance b at the initial time t = 0. 
Then the binding energy is [k n = n~] 
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and the normalised wave function 
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where we note that the physical state is undetermined up 
to a phase (j> which will be employed to find a solution 
of the time dependent Schroedinger equation. We now 
solve the time dependent Schroedinger equation with the 
boundary condition that the wall is moving inwards with 
the speed v. The mathematical problem is 
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together with the initial state (g). Of course the moving 
wall can be equivalently formulated as a time dependent 
5- or Hill- Wheeler potential. This case can be considered 
as a special case of EOI]. The normalised solution reads 
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which determines the phase <j>, especially it leads to the 
initial phase at t = 
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Insertion into (Q) verifies the solution. In |l0| such classes 
of solutions have been used to expand any initial condi- 
tion at t = 0. Here we want to point out that already the 
basic solution (|j) with the extra phase <f> at inital time 
t = bears a physical meaning in itself. The solution 
(f|) represents a non-separable solution. Other classes of 
potentials which admit a separable solution can be found 
in Q. 

The probability density p and current s are easily com- 
puted 
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with p+d x s = as it should. The kinetic energy becomes 
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which shows that even at time t = the initial binding 
energy (@), E° = ^f, is shifted by 
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with the projectile energy E pro] = mv 2 /2 due to the fi- 
nite velocity v corresponding to the finite phase <j). This 
means that due to a finite velocity v of the projectile, 
the target system gets a phase jump <f> and an energy 
immediately when they touch each other. The kinetic 
energy of a level increases with time as 



E n (t) = E^ + E n 
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Ionisation happens when this energy becomes larger or 
equal a now introduced threshold E c such that the ef- 
fective initial binding energy would be E h n = E c - E°. 
Consequently we obtain for ionisation 
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Therefore we find immediate ionisation if the phase ki- 
netic energy E$ exceeds the threshold E c 
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For projectile kinetic energies below we have to wait long 
enough to reach the threshold which leads to 
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where the first inequality comes from the restriction of 
the model in that there should be space between the walls 
left at all. 

Therefore we have two cases. If the projectile is fast 
enough ( |lT| ) , roughly larger than Fermi energy, we knock 
out particles at the first instant of touching. This can be 
seen in analogy to the observation of jjjj where a model 
of instantly removed walls was studied. Due to the time 
dependent solution here we can give the velocity criterion 
where such effect should happen. 

For slower projectiles and fast enough reaction time 
to prevent many - body equilibration and larger reaction 
time than the critical time (12) we will have a knock-out 
of the outer shell nuclei as well. Both cases should be 
considered as a nonequilibrium mode. 

The latter condition for projectiles below or around 
Fermi energy, ([l2]), can be translated into a geometrical 
condition. For the model we assume two approaching 
spherical nuclei with equal radii R. The case of different 
radii is straightforward. The impact parameter B should 
then be smaller than the sum of radii of the two nuclei 
B < 2R in order to allow the necessary overlap. Assum- 
ing constant projectile velocity v pi the time between the 
first touching of the nuclei and the closest approach is 
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Here we approximate the relative velocity v of the nuclei 
by a constant velocity such that at a time t = we have 
the distance R and at the time of maximal overlap we 
have the distance B/2, and the effective velocity v — 
v p (R - B/2)/y/{2R) 2 - B 2 which leads to (|l|). 

The condition for reaction time ( |l2| ) together with 
2R > B translates now into a condition for the geom- 
etry 
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In other words the condition (|l4|) gives the simple re- 
striction on the reaction geometry concerning radii and 
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impact parameter for which an outer bound state char- 
acterised by the binding energy E® is ionised. Rewriting 
( |l4| ) the condition for the projectile velocity reads 
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Therefore we have two projectile energy ranges. Either 
( pi] ) for which ionisation hap pen s immediately or for 
lower projectile energies than ( |lq ) we have overlap and 
possible ionisation if the projectile energy is not too low 
(|l5|). These lower boundary in projectile energy clearly 
distinguishes this model from the diabatic model M. 

In order to draw possible predictions to observable ef- 
fects one would like to have the characteristic angular 
distribution of such emitted particles due to phase space 
squeezing. Since we do not have the exact solution of the 
two dimensional geometry we associate approximatively 
the emission angle with the ratio of impact parameter 
to the distance of nuclei sin a = B/(b — vt) which cor- 
responds to the above geometrical model. The angular 
distribution will be according to the momentum of the 
Wigner function. The latter one can be given with (|J) 
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and ht; — p(b — vt) + mvR. Integrating over all allowed 
spatial coordinates and using the angular relation as de- 
scribed above we obtain the angular distribution which is 
momentum and impact parameter dependent. In figure [l] 
we see that for small impact parameter indeed there is a 
longitudinal emission pattern in agreement with the find- 
ing of the diabatic model |p3|p| . But for specific impact 
parameter around 6 — 8fm for a projectile energy around 
the Fermi energy we see a clear transversal distribution. 
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FIG. 1. The Wigner function ( |l6[ ) integrated over the spa- 
tial coordinates and zero momentum versus angle and impact 
parameter according to b — vt = B / sin (a). The projectile 
energy is around Fermi energy. For nonzero momenta the 
transversal maximum is shifted to lower impact parameter. 

When such ionisation due to phase space squeez- 
ing happens, one can expect small energies beyond the 
threshold. Therefore it should be possible to observe such 
a mode as emitted particles with small total energy and 
the angular distribution should be transversal symmetric 
for the emitted particles. A good plot to observe this 
mode would be the transverse energy versus the ratio of 
the transverse to total energy. The usual multifragmen- 
tation products are on a clear correlation line between 
these two observables. The described mode here should 
be visible as a correlation with small transverse energy 
but a large ratio between transverse and total energy. 

In order to verify the existence of such a mode we 
solved the nonlocal kinetic equation (nonlocal BUU) 
[fill |l6| which leads to figure ||. We have used a soft 
parameterisation of the mean-field and realistic nuclear 
potentials, for details see [0. In the second panel we 
plotted the total kinetic energy of free particles includ- 
ing their Fermi motion. Calculating the current one could 
get rid of the latter motion jlfj] in order to come nearer to 
the experiments which leads of course to a main almost 
linear correlation starting from zero. An alternative way 
would be to use coalescence models which we want not 
to use here in order to maintain theoretical consistency. 
The particles are considered as free if their kinetic energy 
overcomes the binding mean field energy. 

We see a separation of the energy distribution of emit- 
ted particles into two branches at 60fm/c. The main 
line almost at E tr ans/E = 2/3 are the emitted particles 
due to thermal emission and multifragmentation. Besides 
this line we recognise some events in the right lower cor- 
ner of this plot. These 0.2% of emitted particles are due 
to the squeezing mode since they have distinctly smaller 
energy than the rest of emitted particles and are clearly 
transversal. Interestingly this mode is seen at the time 
of closest approach and again later at 120fm/c where the 
neck structure appears. This can be understood since 
during the time in between this mode is shadowed and 
screened by the two nuclei and other emission channels. 

From the right panels one sees that these emitted parti- 
cles originate really from the surface and especially from 
the touching point. The number of indicated trajectory 
points is not representing the amount of emitted particles 
since it is only a cut in the [x, z) plane. The total num- 
ber is given in the middle panel as plot-label. From this 
we see that the homogeneous surface emitted particles at 
the beginning are very small and not due to the squeez- 
ing mode. At the time of remarkable amount of squeezed 
particles one sees that they come mainly from the region 
of overlap. Moreover the potential surface plot shows the 
deepening and squeezing of the potential at the touching 
point. Both observations suggest that these emitted par- 
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tides are probably due to the squeezing mode described 
above. 

Let us comment that the local Boltzmann equation 
(BUU) leads to an even more pronounced effect. Since we 
consider the nonlocal extension of the Boltzmann equa- 
tion as more realistic we give here the smallest estimate of 
the effect. Moreover, the simplified picture given above 
neglects completely the rebinding mechanisms e.g. by 
mean fields which will limit the ionisation. This was a 
sacrifice for analytical solvability. However the realistic 
simulation shows that a small percentage of events might 
show this behaviour. 

To summarise a new mode is predicted due to fast 
shrinking of available phase space for outer particles 
when two nuclei collide. The time in such collisions is 
too short to establish a common ionisation threshold. In- 
stead the energy of the level increases and ionisation can 
occur. In opposition to the Mott transition, here the 
levels increase which gives rise to the name "squeezing 
mode" . The characteristics of such a mode are that the 
emitted particles have very small total energy and are 
transversal. This feature and the fact that the mode ap- 
pears at the beginning of the collision distinguishes it 
from the "towing mode" p8| . In a corresponding plot we 
could identify such events in nonlocal BUU simulations. 
In order to demonstrate this mode an analytical solution 
of the time dependent Schroedinger equation has been 
given. 

The described mode is not restricted to nuclear colli- 
sions. Anytime when two clusters of particles collide fast 
enough to prevent the formation of a common threshold 
and long enough to squeeze the outer bound states this 
squeezing mode should be possible to observe. 

The discussions with Jean Colin, Daniel Cussol and 
Jaques Normand are gratefully acknowledged which have 
convinced us that this mode is possible to observe experi- 
mentally. Especially to Victor Dodonov, Marco Ameduri 
and John Frankland we are indebted for valuable com- 
ments. 
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FIG. 2. The evolution of a Ta + Au collision at 33 MeV lab energy and 8 fm impact parameter. The left panel gives the 
spatial density contour-plot in the (x, 0, z) plane and the local currents as arrows. The middle panel shows the transversal 
energy E t distribution versus the ratio of transversal to total energy Et/E. The amount of emitted particles due to the squeezed 
mode which is located in the right lower corner is given as plot label. The right panel represents the contour-plots of the total 
energy of particles. The dark areas gives the positive total energy indicating areas where particle can escape and the lighter 
colour scales the deepness of negative total energy indicating bound states. The black dots are the position of test particles in 
the (x, 0, z) plane which contribute to the new squeezing mode. 
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